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Essence of inviscid shear instability: a point view of vortex dynamics 
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The essence of shear instability is fully revealed both mathematically and physically. A general 
sufficient and necessary stable criterion is obtained analytically within linear context, ft is the 
analogue of Kelvin- Arnol'd theorem, i.e., the stable flow minimizes the kinetic energy associated 
with vorticity. Then the mechanism of shear instability is explored by combining the mechanisms 
of both Kelvin-Helmholtz instability (K-H instability) and resonance of waves. It requires both 
concentrated vortex and resonant waves for the instability. The waves, which have same phase 
speed with the concentrated vortex, have interactions with the vortex to trigger the instability. 
We call this mechanism as "concentrated vortex instability". The physical explanation of shear 
instability is also sketched. Finally, some useful criteria are derived from the theorem. These results 
would intrigue future works to investigate the other hydrodynamic instabilities. 

PACS numbers: 47.15.Fe, 47.20.-k, 47.32.-y 



The hydrodynamic instability is a fundamental prob- 
lem in many fields, such as fluid dynamics, astro dynam- 
ics, oceanography, meteorology, etc. There are many 
kinds of hydrodynamic instabilities, e.g., shear instability 
due to velocity shear, thermal instability due to heating, 
viscous instability due to viscosity and centrifugal insta- 
bility due to rotation, etc. Among them, the shear insta- 
bility is the most important and the simplest one, which 
has been intensively explored (see Q, 0, Q and references 
therein). Both linear and nonlinear stabilities of shear 
flow have been considered, and some important conclu- 
sions have also be obtained from the investigations. 

On the one hand, the nonlinear stability of shear How 
has been investigated via variational principles. Kelvin 
and Arnol'd H, IE have developed variational prin- 
ciples for two-dimensional inviscid flow § . They showed 
that the steady flows are the stationary solutions of the 
energy H . And if the second variation S 2 H is definite, 
then the steady flow is nonlincarly stable. Moreover, 
Arnol'd proved that the flow is linearly stable provided 
that 5 2 H is positive definite, and he also proved two non- 
linear stability criteria [a, S H 01 ■ However, d 2 H is al- 
ways indefinite in sign, except for two special cases (see 
and references therein). Though variational principle 
is indeed powerful, it is also inconvenient for real appli- 
cations due to the lack of the explicit expressions in both 
H and the stability criteria. 

On the other hand, the linear stability of shear flow has 
also been investigated via Rayleigh's equation. Within 
the linear context, there are three important general 
stability criteria, which are Rayleigh's criterion [To| . 
Fj0rtoft's criterion 01 an d Sun's criterion 0. As all 
the criteria have explicit expressions, they are more con- 
venient in real applications, and are widely used in many 
fields. Based on the former investigations, Sun 0] also 



pointed out that the flow is stable for Rayleigh's quotient 
/(/) > (see Eq.Q behind). However, the sufficient cri- 
terion for instability still lacks. 

To understand the shear instability, some mechanisms 
were suggested. Among them, Kelvin-Helmholtz insta- 
bility (K-H instability) is always taken as a prototype, 
which isphysically explained as the instability of a sheet 
vortex |1J| . An another mechanism of instability is due to 
the resonance of waves 0, El El El- Butler and Far- 
rell |l5j| clearly showed with numerical simulations that 
the resonance introduces an algebraic growth term into 
the temporal development of a disturbance. Baines and 
Mitsudera 0] also used broken-line profile velocity as 
a prototype to explain the interaction of waves. Their 
explanation is so brilliant that it can explain why the 
instability occurs for a finite range of wavenumber and 
how the waves amplify each other. However, both mecha- 
nisms are independent of the basic flows, for that Kelvin- 
Helmholtz model deals only with vortex and resonance 
mechanism only considers the waves. Thus the relation- 
ships between those mechanisms and the shear flows are 
still covered. 

Overview the former investigations, the essence of the 
shear instability is remained to be elucidated. Though 
H is very important in variational principles, both the 
explicit expression and the meaning of it have not been 
revealed before. The connection between linear and non- 
linear stability criteria should be retrieved explicitly. The 
relationships between instability mechanisms and the 
physical explanation for shear instability are needed. The 
aim of this letter is to fully reveal the essence of the shear 
instability by investigating the inviscid shear flows in a 
channel. And other instabilities in hydrodynamics may 
also be understood via the investigation here. 

For the two-dimensional inviscid flows with the velocity 
of U, the vorticity £ = V x U is conserved along pathlines 
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Its linear disturbance reduces to Rayleigh's equation pro- 
vided the basic flow U being parallel. Consider an shear 
flow with parallel horizontal velocity U(y) in a channel, as 
shown in Fig^ The amplitude of disturbed flow stream- 
function ip, namely 4>, satisfies 0,0,13 : 



k 



U" 



u- 



= 0, 



(2) 



where k is the nonnegative real wavenumber and c = 
c r + ici is the complex phase speed and double prime 
" denotes the second derivative with respect to y. The 
real part c r is the phase speed of wave, and Cj ^ de- 
notes instability. This equation is to be solved subject to 
homogeneous boundary conditions 4> = at y = a, b. 

Consider that the velocity profile U (y) has an inflection 
point y s at which U" = U"{y s ) — and U s — U{y s ). 
As Sun 01 nas pointed out, the following Rayleigh's 
quotient 1(f) > implies the flow is stable. 



1(f) 



I b a [H'\\ 2 + f(y)H\\ 2 }dy 
JaUW 2 



(3) 



where f(y) 



u" 
u-u„ 



While if /(/) < 0, then there is 



a neutral stable model with k 2 ^ = —1(f) and c r = U s . 
Moreover, there are unstable modes with c r = U s and 
Cj 7^ if /(/) < and < k < kpf. This can also be 
proved by following the way by Tollmien |l8j . Friedrichs 

llHlSl and Lin El- Thus !(f) = means the now 
is neutrally stable. And there is only one neutral mode 
with k = and c r = U s in the How. These conclusions 
can be summarized as a new theorem. 

Tollmien-Fridrichs-Lin theorem: The flows are neu- 
trally stable, if /(/) = 0. The flows are stable and un- 
stable for /(/) > and /(/) < 0, respectively. 

From the above proof, it is obvious that the shortwaves 
(e.g. k ^> kjy) are always more stable than longwaves 
(e.g. k <C fcjy) in the inviscid shear flows USUI. So the 
shear instability is due to long-wave instability, and the 
disturbances of shortwaves can be damped by the shear 
itself without any viscosity |22^ . 

As mentioned above, the linear stability criterion can 
be derived from the nonlinear one 6|. Moreover, the 
nonlinear stability criterion can also be obtained from 
the linear one. To illuminate this, the nonlinear criterion 
is retrieved explicitly via the above theorem, which is 
briefly proved as follows. 

Similar to Arnol'd's definition, the general energy H 
here is defined as 



H 



.1, 



h(9)]dy 



(4) 



where ^(y) is the streamfunction of the flow with ^ = 
U (y) — U s , and h(^>) is a function of 'J. The variation of 
5H = gives 



A* = h'(^). 



(5) 



So Arnol'd's variational principle is retrieved. The func- 
tion h = ||V*|| 2 /2 can also be solved from Eq.JSJ). The 
second variation S 2 H holds 



S 2 H 



[\\y^\\ 2 + h"(^}dy, 



(6) 



where ip denotes the variation 8^. Noting that h"(^) 
remains unknown, Arnol'd's nonlinear criteria have not 
be extensively used. Fortunately, we can obtain the ex- 
plicit expression of h"{^) here via the investigation on 

<>y 



the linear stability criterion. For IP =U — U s , 



, »/ t s dh' db! dy 



U" 



dy cM (U - Us 



= M- (7) 



As h"(^>) is solved explicitly, 8 2 H has an explicit expres- 
sion, which is greatly helpful for real applications. 

Let the streamfunction of the perturbation ip(x, y, t) = 
4>(y)e l ( kx ~ ut \ where u is the frequency. The averages of 
||VV>|| 2 and HV-II 2 are (||<//|| 2 +fc 2 ||</)|| 2 )/2 and ||0|| 2 /2 along 
the How direction x, respectively. So Eq.© reduces to 



S 2 H 



f(y)U\\ 2 ]dy. (8) 



The sign of 8 2 H is then associated with /(/) in Eq.Q. If 
1(f) < 0, the second variation 8 2 H can be both negative 
and positive, i.e., the stationary solution is a saddle point. 
And /(/) > implies 8 2 H is positive definite and vice 
versa. So the stable flow has the minimum value of the 
total kinetic energy H . The physical meaning of H can 
also be revealed, as the explicit expressions of H , ^ = 
U - U s and h"(^) = f(y) are obtained. 

First, according to the expressions, the velocity U in 
vorticity conservation law Eq.Q can be decomposed to 
two parts: the rotational flow U — U s and the irrotational 
advection flow U s . The vorticity £ in Eq.QJ depends 
only on U — U g , and U s is only advection velocity. Then 
U — U s and U s are associated with the dynamics and 
kinetics of the flow, respectively. Eq. physically shows 
that the vorticity field is advected by U s , which can also 
known from the conservation of vorticity in the inviscid 
flows. A similar example is the dynamics of vortex in the 
wake behind cylinder, where the vortices dominate the 
dynamics of the flow and they are advected by mean How 
(see Fig. 2 in |23). The decomposition of velocity may 
be useful in vortex dynamics, for that our investigation 
clearly shows that the dynamics of the flow is dominated 
by vorticity distribution. 

Then the physical meaning of H can also be under- 
stood from the above investigation. It is not U but 
U — U s that is associated with the general energy H, so 
H = I f b (U — U s ) 2 dy is not the total kinetic energy but 
the kinetic energy of flow with vorticity. Thus the stable 
steady states are always minimizing the kinetic energy of 
the flow associated with vorticity. This is also the reason 
why the flow with maximum vorticity might be unstable, 
as Fj0rtoft's criterion shows. We would like to restate it 
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as a theorem and to name it after Kelvin jj| and Arnol'd 
6J for their contributions on this field @. 

Kelvin- Arnol'd theorem: the stable flow minimizes the 
kinetic energy of flow associated with vorticity. 

Both Tollmien-Fridrichs-Lin theorem and Kelvin- 
Arnol'd theorem are equivalent to the following simple 
principle [24|]: The flow is stable, if and only if all the 
disturbances with c r = U s are neutrally stable. 

We have obtained the sufficient and necessary condi- 
tions for instability, then the physical mechanism of in- 
stability can be understood from them. The following 
investigation will reveal that the essence of shear insta- 
bility is due to the interaction between the " concentrated 
vortex" and the corresponding resonant waves. 

The concept of "concentrated vortex" comes from 
Fj0rtoft's and Sun's criteria, where the necessary condi- 
tions for instability require that the base vorticity must 
be concentrated enough 0. We call it "concentrated 
vortex" for latter convenience. The concentrated vor- 
tex is a general model of sheet vortex in the Kelvin- 
Helmholtz model, for that the sheet vortex can be recov- 
ered as the concentrated vortex £(y s ) — * oo (see [2^ for 
a comprehensive discussion about the Kelvin- Helmholtz 
model and continued shear profiles). Then how the shear 
flow becomes unstable, if there is a concentrated vortex? 
As the sufficient condition for instability is 1(f) < 0, 
the normal modes in the regime of < k < with 
Cr = U s are unstable. They are stationary or standing 
waves, comparing to the velocity at inflection point. So 
that the shear instability is due to the disturbance of con- 
centrated vortex by the standing waves with c r — U s . In 
this case, the resonance mechanism is valid. The inter- 
action waves propagate at the same speed with the con- 
centrated vortex, so that they are locked together and 
amplified by simple advection |16|. In short, the distur- 
bances on concentrated vortex is amplified like that in 
Kelvin-Helmholtz model. 

This instability mechanism combines both K-H insta- 
bility and resonance mechanism. Physically, the standing 
waves (with c r = U s ) can interact with the concentrated 
vortex, so they can trigger instability in the flows. While 
the travelling waves (with c r U s ) have no interaction 
with the concentrated vortex, so that they can not trigger 
instability in the flows. This is the mechanism of shear 
instability. As pointed out above, the shear instability is 
due to long-wave instability. If the longwaves are unsta- 
ble, they can obtain the energy from background flows. 
In this way, the energy within small scales transfer to 
and concentrate on large scales. In this way, the shear 
instability itself provides a mechanism to inverse energy 
cascade and to maintain the large structures or coherent 
structures in the complex flows. 

To illuminate the mechanism of shear instability, a 
physical diagram is also presented here by following 
the way of interpreting the K-H instability Figffl 
sketches the mechanism of shear instability in terms of 
wave disturbances of vortices. The mean velocity pro- 
file U(y) has an inflection point at y s = with U s = 0, 
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FIG. 1: Sketch of shear instability: physical interpretation. 
Left parts depict the profiles of velocity U(y) and vorticity 
£(y), right ones depict the disturbance of vorticities. The un- 
stable veloctiy profile U (y) has a local maximum in vorticity 
£(y) (a). If the vortices (A, B and C) disturbed from their 
original positions (dashed line) to new places (solid curve), 
they will be taken away from their original positions due to 
pressure difference. The stable veloctiy profile U(y) has a lo- 
cal minimum in vorticity £(j/) (b). The disturbed vortices (A', 
B' and C) will be brought back to their original positions due 
to pressure difference. 

and the corresponding vorticity is t;(y) = —U'(y). There 
is a local maximum at y s in the unstable vorticity pro- 
file in Fig^,. According to Eq.JTJ, the vorticity is con- 
served in the inviscid flows. If the vortices at the local 
maximum (A, B and C) are sinusoidally disturbed from 
their original positions (dashed line) to new places (solid 
curve), they have negative vorticities with respect to the 
undisturbed ones. The vortices will induce cyclone flows 
around them in consequence. The flows around the vor- 
tices become faster (slower) in the upper (lower) of vortex 
A referred to the basic flow U(y). The pressures at up- 
per and lower decrease (indicated by - signs) and increase 
(indicated by + signs) according to Bernoulli's theorem, 
respectively. Then vortex A gets a upward acceleration 
due to the disturbed pressure difference as the uparrow 
shows. This tends to take the vortex away from its orig- 
inal position, so the flow is unstable. On the other hand, 
FigJIJ) depicts the disturbances in a stable velocity pro- 
file, where a local minimum is in the vorticity profile. 
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The disturbed vortices have positive vorticities with re- 
spect to the undisturbed ones. The vortices will induce 
anticyclone flows around them in consequence. So vor- 
tex A' get a downward acceleration due to the disturbed 
pressure difference as the downarrow shows. This tends 
to bring the vortex back from its original position, so 
the flow is stable. In this interpretation, the advection 
of U s is independent of the shear instability, only the 
flow field of U — U s and the corresponding vorticity £ 
are the dominations. The unstable disturbances in Fig 
n have c r — U s , which consists with the former discus- 
sions. This physically explains why the maximum and 
minimum vorticities have different stable aspects. 

Though Tollmien-Fridrichs-Lin theorem is a sufficient 
and necessary condition for stability, the unknown <fi in 
the theorem restricts its application. So some simple cri- 
teria may be more useful. For the parallel flows within 
interval a < y < b, there are two simple criteria. 

Corollary 1: The flow is stable for f(y) > -(^) 2 
[30. 
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